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ABSTRACT 

We present a detailed analysis of an astrophysical mechanism that generates cosmo¬ 
logical magnetic fields during the Epoch of Reionization. It is based on the photoion¬ 
ization of the Intergalactic Medium by the first sources formed in the Universe. First 
the induction equation is derived, then the characteristic length and time scales of the 
mechanism are identified, and finally numerical applications are carried out for first 
stars, primordial galaxies and distant powerful quasars. In these simple examples, the 
strength of the generated magnetic fields varies between the order of 10“^^ G on hun¬ 
dreds of kiloparsecs to 10“^® G on hundreds of parsecs in the neutral Intergalactic 
Medium between the Stromgren spheres of the sources. Thus this mechanism con¬ 
tributes to the premagnetization of the whole Universe before large scale structures 
are in place. It operates with any ionizing source, at any time during the Epoch of 
Reionization. Finally, the generated fields possess a characteristic spatial configuration 
which may help discriminate these seeds from those produced by different mechanisms. 

Key words: Gosmologyrtheory - dark ages, reionization, first stars - large-scale 
structure of Universe - magnetic fields - intergalactic medium. 


1 INTRODUCTION 

The origin of cosmological magnetic fields is a major open 
question in cosmology. Recent high energy gamma ray obser¬ 
vations suggest that a substantial fraction, if not the whole, 


of the intergalactic space is magnetized (e.g. Neronov & 
|Vovk 2010|. The current paradigm to account for the ex¬ 


istence of such cosmological magnetic fields states that a 
first mechanism, or several mechanisms combined, gener¬ 
ated large scale magnetic fields but of very weak strengths 
(so-called seed fields) that were amplified later on, during 
structure formation, essentially through turbulence (see, for 
instance, [Brandenburg fc Subramanian||2005| |Ryu et al.| 


2011 Widrow et al.|20TT and references therein). Numerous 


mechanisms for generating magnetic fields on cosmological 
scales have been proposed, operating mainly in the primor¬ 
dial Universe, during inflation or the electroweak and quark- 
hadron phase transitions (e.g. Grasso & Rubinstein 2000[ 


Widrow 2002 Kandus et al. 2011 and references therein). 


However, mechanisms operating in the radiation dominated 


era (Harrison 1970 

Zakharov & Anikanov 

19921 or during 

recombination (e.g. 

Berezhiani & Dolgov 

2004| Takahashi 

et al.|2008| Fenu et al. 20111, requiring some level of vorticity 


(possibly re-generated at the second order in perturbations). 


E-mail:jdurrive@ias.u-psud.fr; 


have also been proposed. Finally, astrophysical processes op¬ 
erating after recombination, capable of generating magnetic 
fields of cosmological interest have been investigated too. 
They include the well-known Biermann battery, due to a 
thermal electromotive force, that was first introduced in the 
context of stars ( Biermann||1950 1 and later successfully con¬ 
sidered in cosmological contexts such as structure formation 


i Pudritz fc Silk||l989 Kulsrud et al. jl997 Xu et al. 20081 


and the propagation of ionization fronts | Subramanian et al. 
1994 Gnedin et al.|2060 1 during Cosmological Reionization. 
Gollision-less shocks in cosmology are also potentially capa¬ 
ble of generating magnetic fields by triggering plasma in¬ 
stabilities (e.g. Medvedev et al.||2006 Coroniti||2014 1. Note 
that large scale magnetic fields may actually have been cre¬ 
ated within galaxies and then transported into the inter¬ 


galactic medium by powerful winds and/or jets (Kronberg 
et al.[1999 Furlanetto fc Loeb|2001 Beck et al.|2013 1. Glob¬ 
ally, the level at which all these mechanisms may have con¬ 
tributed to the magnetization of the intergalactic medium 
is still an open question (for recent reviews, see for example 


Widrow et al.|[^011[ [Durrer fc Neronov||2013l |Subramanian| 


20151 to which future observations with radio interferome¬ 


ters will bring essential pieces of answer (e.g. Bonafede et al. 


Revisiting [Langer et al.| ( |2005[ ) , we present here a de¬ 
tailed analysis of an astrophysical mechanism, based on the 
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photoionization of the Intergalactic Medium (IGM), bound 
to have operated during the first billion years of the Uni¬ 
verse. Ando et al. (20101 and Doi & Susa (2011) explored 


numerically the same mechanism and compared it to the 
Biermann battery. However, in their analysis, they focused 
on the competition between these two mechanisms at the 
boundaries of self-shielded, essentially neutral clumps em¬ 
bedded inside the Hll regions of individual first stars. In 
these conditions, they found that the Biermann battery pro¬ 
duces stronger magnetic fields than the radiation effects, on 
hundreds of parsecs scales. Here, we study analytically this 
mechanism, relying on the momentum transfer from ioniz¬ 
ing photons to electrons, on large scales, way outside the 
Stromgren spheres of clusters of Population HI stars, pri¬ 
mordial galaxies and quasars. 

In this paper, in sectionwe introduce the mechanism 
in full generality. Then, in section]^ we obtain a simplified 
order of magnitude estimate of the magnetic field strength 
as well as convenient, although approximate, scaling rela¬ 
tions. In section]^ we analyse in full details the expression 
obtained in section apply it to the context of the Epoch 
of Reionization (EoR) and obtain numerical values of the 
magnetic fields generated in the IGM by the first luminous 
sources. Finally, section|^is dedicated to a discussion where 
a comparison with the Biermann battery is also included. 


2 THE MECHANISM 
2.1 Presentation 


2.2 Formalism 

2.2.1 Fields 


In the non relativistic limit, displacement currents are neg¬ 
ligible, so we consider the following Maxwell’s equations: 

V X 5 = —J V-.B = 0 

c 

^xE = --dtB V-E = p 

c 

where J is the total current density and p the total charge 
density. We will take p = G since the characteristic length 
scales of the problem are much larger than the Debye length, 
of the order of the kilometre here. Initially the current, the 
electric field and the magnetic field are null since we are 
interested in their ab initio generation. 

Cosmological recombination was an incomplete process: 
during the Dark Ages, a tiny non zero ionization fraction 
remained in the IGM. We will call residual electrons and 
residual ions the free electrons and ions from this plasma. As 
sources switched on, they photoionized their surroundings 
thus liberating new electrons and ions. As we will see, it is 
instructive to consider separately these two types of charged 
species. Hence we will consider five different species, namely: 


a = < 


1 : residual electrons; 

2 : residual ions; 

3 : photoionization electrons; 

4 : photoionization ions; 

5 : neutrals. 


The first sources in the Universe switched on in an essen¬ 
tially neutral IGM, mostly made of Hydrogen, below red- 
shifts of 30 (e.g. Loeb fc Furlanetto|2013 1. As they radiated, 
the sources formed fully ionized regions around them, called 
Stromgren spheres, created mainly by photons just above 
the ionization threshold of 13.6 eV. However, higher en¬ 
ergy photons were able to escape the Stromgren spheres and 
propagate deeper into the IGM, because the photoionization 
mean free path is proportional to the cube of their energy. 
As pointed out by |Langer et al.| ( |2005[ ), these photons trans¬ 
ferred their momentum to electrons in the surrounding, oth¬ 
erwise neutral medium, and thus generated radial currents. 
These currents were in turn able to induce large scale mag¬ 
netic fields, provided the corresponding electric fields were 
rotational. This condition was actually satisfied thanks to 
the anisotropic absorption of the radiation due to the inho¬ 
mogeneities of the neutral IGM. 

Formally speaking, the ionization process is described 
microscopically as a perturbation of the distribution func¬ 
tion of electrons. The description is then reduced to a macro¬ 
scopic monofluid description to get a generalized Ohm’s law, 
the rotational of which leads to the induction equation (sec¬ 
tion 2.2). The general expression for the generated magnetic 
field thus obtained is then examined in a simple model of 
the cosmological context we are interested in, using power 
law spectra for the sources and Gaussian profiles to model 
the dumpiness of the IGM. This allows us to identify the 
characteristic properties (characteristic length scales, typi¬ 
cal strengths generated and field lines) of the regions that 
are significantly magnetized (section 4.3), and then to ob¬ 
tain numerical estimates of these photogenerated magnetic 
fields (section |4.4[). 


Each of these matter fields is characterized by its distri¬ 
bution function fa. Each fa is governed by the following 
generalized Vlasov equation: 


or I ^ ^fct Qa 

dtfa+V.^;— -I- 

or rria 


E + 


V X B 


dfa 

dv 


dtfa\^+ dtfa\^ , 


( 1 ) 

where qa is the charge of species a and rria its mass. On the 
right hand side, the first term is the usual collision term and 
the second is the source term due to photoionizations that 
is detailed in appendix]^ 

Astrophysical ionizing sources are characterized by their 
specific spectral density E. In principle one should solve the 
complete radiative transfer equation governing the evolu¬ 
tion of E, but for our purpose it is enough to consider the 
following solution to this equation: 


E{t,r,k) = L. ^S{k-f), (2) 

where is the spectral luminosity density of the source 
and Tu = (Ti^ riHidr is the optical depth with a^, the pho¬ 
toionization cross section. The S(k — f) factor accounts for 
the fact that the source emits radially and corresponds 
to the geometric dilution. 


2.2.2 Induction equation 

This kinetic description is convenient to build the right- 
hand-side of equation ([^, but it contains much more in¬ 
formation than needed since we are looking for the macro¬ 
scopic and summed over all species quantity J appearing 
in Maxwell’s equations. Thus we reduce our description to 
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that of a monofluid (cf. Krall & Trivelpiece 1973 and ap¬ 


pendix and get the following equation governing J, the 
generalized Ohm’s law: 


dtJ+ (V.V'J J+(v-jjV-VV-Vp 




-P + C+Eo 


where 


Pa = y mavdlfa\,d^V 


-Pc, 


( 3 ) 


( 4 ) 


and the other terms are detailed in appendix The last 
term of equation corresponds to the momentum transfer 
from photons to electrons. 

In this paper, we are interested in the generation of mag¬ 
netic fields in a cosmological context. The typical gradient 
scales L correspond to the size of matter inhomogeneities in 
the high redshift intergalactic medium. To get an idea, we 
can consider L ~ 10 kpc. Further, as shown in |Langer et al.| 
(20051, on the very short initial times typical of the plasma 


time-scales, the generated strengths of the fields are negli¬ 
gible, and the next characteristic time is set by the lifetime 
of the source, typically from 1 to 100 Myrs during Reioniza¬ 
tion. Therefore, when needed, we will take T ~ 10 Myrs as 
a typical time-scale. Finally, during Reionization the typical 
residual electron density is rie ~ 2 x 10“^h = 4 x 10“® cm“ 


at 2 = 9 (e.g. Seager et al.||20(io| . We will use these values 
in the following order of magnitude estimates. 

The general expression (|^ may then be simplified as 
follows. 

- The time variation of J is completely negligible here. 
Indeed, combining Maxwell’s equations yields - 
that: 


c'^T 


-E 


dtJ 


2 

o 

meC^ 


10^® > 1 


( 5 ) 


where q is the charge of an electron. 

- Since displacement currents are neglected, J has a ro¬ 
tational form from Maxwell’s equations, so V • J = 0 here. 

- Starting from zero, J, E, and B are first order terms 
initially and we can therefore linearize equation §. Note 
that the generated magnetic held strengths will be small 
enough for this linearization to remain valid during the time- 
scales of interest. 

- Fluids with protons are assumed to move slowly. Thus 
Va for ionic species are first order terms. Since J and B are 
also of first order and me ^ rrii , the second term on the right 
hand side of ([^ is a second order term and is neglected. 

- For the pressure term, we assume the V ■ for pro¬ 
tonic species (a = 2,4, and 5) is small due to their large 
inertia. For residual electrons, we neglect viscosity effects en¬ 
suring an isotropic pressure tensor. Photoelectrons, however, 
are generated radially thus introducing a small anisotropy, 
but we neglect this with respect to the pressure gradients of 
residual electrons. Hence: 

P=-—Vpe (6) 

rUe 

where Pe is the residual electron pressure. 

- For the collision term, numerous types of collisions 
could in principle be considered given all the species in¬ 
volved. However, it would prove unnecessary since all the 


collision frequencies are far too small and therefore the as¬ 
sociated time-scales far too large with respect to T. More 
precisely, taking the usual linear approximation C = VcJ 
where Uc is an averaged collision frequency, and comparing 
this term for example to the electric field in the plasma yields 


E 


rrie 


c 


q^rie D vp 
rUe T 




Since typical values for collision frequencies hardly exceed 
Hz, this ratio is always huge. 

- In the last term of equation ([^, the dominant and es¬ 
sential contribution comes from the momentum exchange 
between ionizing photons and photoelectrons. Other contri¬ 
butions are negligible because of the large inertia of pro¬ 
tons and neutrals, and the Thomson scattering cross section 
which is many orders of magnitude smaller than the pho¬ 
toionization cross section. This last term therefore reduces 
to 


/ 


Ppe = / rUeV dtfpel^av 


( 8 ) 


where ‘pe’ stands for photoelectrons. This expression may 
be interpreted as follows: dtfpe\g d^v d^f dt is the number of 
photoelectrons generated in a volume element during dt, 
appearing with momentum rrieV. Thus Ppe d^f dt represents 
the total electron momentum appearing in a volume d^f 
during dt. Equation (1^ is a momentum density generation 
rate. While equation has been correctly described many 
times as a close analogue to Newton’s second law, we stress, 
however, that the term ([^ is not, in essence, a force density, 
but a source of momentum. Further, as detailed in appendix 
[A{ expression (|^ may be explicited as 



where /mt is the fraction of momentum transferred from a 
photon to an electron in the photoionization process, and zzq 
is the Hydrogen ionization threshold. 


Finally, under these assumptions. Ohm’s law sim¬ 
plifies to: 

0 =—qn^E — Vpe + Ppe- ( 10 ) 


Note that rie is the total number density of electrons, but 
the source term ppe is only due to the newly photoionized 
electrons, and the pressure term is only due to residual elec¬ 
trons. Furthermore, we emphasize again that this equation 
does not correspond to the balance of forces acting on sin¬ 
gle electrons | |Ando et al.||20fo| |Doi fc Susa||20lT| , but it 
rather expresses the readjustment of the electric field in the 
plasma in response to the apparition of new currents from 
photoelectrons. 


The induction equation is then given by the curl of (101, 
and may be written, using Faraday’s law, as 


= X Vpe 


+ 5 ^ X + VXT nmdr x 4,e) 


( 11 ) 


where Xe = is the total electron fraction, and qpe has the 

same expression as ppe in 0 but with instead of cr^ in the 
integrand. The first term on the right hand side is the usual 
Biermann battery term and the two additional terms are due 
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to photoionization. The Biermann term will be discussed in 
section and will not be considered here otherwise. Then 
integrating the magnetic held at time t and position r 
may be written as a sum of two contributions: 


— Blocal “t“ Bglobal 

where the ‘local’ term is 

Bioca,= f dt 

Jo 


( 12 ) 


(13) 


11/'°° 

Xo /mta.L.e (14) 


qx; 

= VXe X f 

and the ‘global’ term 

^global = / ^ global C’gioijal 

Jo 

1 r / „ 2 T -T, 


Bg,„hal = / K"ob.,F5,!°“, dt 

F 


(15) 


(16) 


- global A o 

qXe 47rr^ 


1 

/ fuit<xtL^e (17) 

J l/Q 


^geom _ y 
-^global ~ ^ 




X r. 


(18) 


Formally, Biocai and Bgiobai are both products of two terms, 
integrated over time: an ‘interaction’ term and a ‘geo¬ 
metric’ term The interaction term characterizes the 

impact of the source at a time t and a position r as it includes 
the absorption, geometric dilution, the photoionization cross 
section and the fraction of momentum transferred from pho¬ 
tons to electrons. The geometric term determines whether 
gradients in the IGM are indeed non radial as required, in¬ 
dependently of the properties of the ionizing source. 


Physically, one can interpret (121 in the following way. 


The charge separation induced by photoionizations gener¬ 
ates an electric held in the plasma satisfying the equilibrium 
(101. For a magnetic held to grow out of it, the electric held 
must have a non vanishing curl. 

To see how this condition may be fulhlled, consider 
two adjacent volume elements at a given distance from the 


source. In both volume elements, the equilibrium (101 is sat- 


ished and dictates the value of the electric held there. For¬ 
getting about the pressure term in this discussion, we see 
that E depends on two things: the local density of electrons 
Tie, and the local ability Ppe of the source to photoionize the 
medium. Thus, there are three ways for E to be different 
in the two volume elements: (i) they have the same Ppe but 
diherent Ue, (ii) they have the same Ue but different ppe, 
and (iii) both rie and Ppe are diherent. Situation (i) may 
occur if the electron density of the plasma is locally inho¬ 
mogeneous. The resulting magnetic held is the local term 
(13l. Situation (ii) occurs when the incident distribution of 
photons is inhomogeneous, that is when the radiation held 
is anisotropic. But since the source itself emits isotropically 
by assumption, the only origin of such anisotropies are inho¬ 
mogeneities in the medium in which photons propagate, so 
that absorptions along adjacent lines of sight diher. There¬ 
fore, (ii) occurs due to non radial gradients of the column 


densities, which gives rise to the global term (161. Finally, 


situation (iii) corresponds to the general case in which every¬ 
thing is inhomogeneous and yields the total magnetic held 

(^. 

Note that if at some distance two adjacent lines of sight 
diher, they will in general remain diherent further away from 


the source. For this reason the global term generates mag¬ 
netic helds on large distances. Therefore, behind an inho¬ 
mogeneity, some magnetic held is generated from this global 
term even if the medium is homogeneous there, and is only 
attenuated by geometric hux dilution, absorption and the 
1/r factor from the gradient. 


3 ORDER OF MAGNITUDE 

First, to get an intuition of the efhciency of this mecha¬ 
nism, let us provide a crude estimate of the reached mag¬ 
netic strengths by evaluating expression (161 in the typical 


case of a primordial galaxy at 2 = 10. For this, as detailed in 
sections |4.1| and |4.2[ we consider the following parameters. 
We assume that the source has a Stromgren radius ~ 48 
kpc, a spectral index a = —2, a lifetime ta = 100 Myrs, and 
a total UV luminosity LtoY = L,ydu = 2.5 x lO^L©. At 
those epochs, the ionized fraction of the IGM is Xe — 2x 10“"* 
while its mean density is h ~ 2.5 x 10“"* cm“®. For the gen¬ 
eration of electric helds possessing a curl, we consider that 
an inhomogeneity of amplitude 5o = 4 is positioned at a dis¬ 
tance D ~ 1.3rs from the source. For illustration, the mag¬ 
netic held strength is evaluated at a distance r = Xs + Iovq 
from the source, where = 27 is the mean free path 
of photons of energy three times the Hydrogen ionization 
threshold. Finally, the strength of the magnetic seed keeps 
on growing linearly in time as long as the momentum ex¬ 
change process between photons and initially bound elec¬ 
trons is active. The photons we are interested in here being 
way above the Hydrogen ionization threshold, they do not 
induce any signihcant expansion of the Stromgren sphere 
(see discussions on tree and ti in the next section), but in¬ 
stead propagate outside deep into the neutral IGM. Thus, 
the limiting time-scale for the process to operate is essen¬ 
tially the life-time of the ionizing source ta- 


To get a simple expression of (161, let us now take ad¬ 


vantage of the strong dependence on v in the optical depth to 
approximate absorption in the interaction term by a Heavi¬ 
side function: ~ 6[v—Vc), where Vc is such that = 1. 

Then 


Va{t, r) = Vo 


/ riHidr^ 


1/3 


~ Vo 


1/3 


(19) 


where the approximation assumes a homogeneous back¬ 
ground. Gonsidering in addition a power law spectrum as 


in (221 on [z/q, oo[ and using the fact that r — ^ 1^,^ yields 


- global ' 


8 1 
35 qXa ^ 


- 7/3 


( 20 ) 


Finally, to estimate the contribution of the geometrical 
term, we consider a spherical inhomogeneity of characteristic 
size a and amplitude So - Then |Egi°b™il = ^dg f nmdr, where 
only the inhomogeneous part nSo of the density matters, 
the angular variation is of the order of the angular diameter 
a/D, and the integral along the line of sight is essentially 
the size of the inhomogeneity a. Hence IF’giob™!! ~ nSo^- We 
may then plug in the numerical values mentioned above and 
get 


|B| ~ 8xl0~^° G ( ^ 


2.5 X lO^L© 


ta 


10® yrs 


• ( 21 ) 
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This rapid estimate provides interesting information on the 
typical magnetic strength one might expect from this mech¬ 
anism. In particular, it shows that momentum transfer from 
photons to electrons during Reionization is capable of pro¬ 
ducing fields with similar intensities as the usual ther¬ 
mal Biermann battery. However, the strength of the held 
strongly depends on the geometry and position of the inho¬ 
mogeneity. Angularly, it is essentially governed by the geo¬ 
metrical term and radially by the interaction term, so that 
a strength of this order is reached only in specihc regions. 

Additionally, in accord with intuition, expression (211 
also shows that stronger gradients (i.e. larger inhomo¬ 
geneities) will create stronger helds, and that longer-lived 
sources are better at creating stronger magnetic seed helds. 
Moreover, it indicates that the reached magnetic strength 
grows linearly with the UV luminosity of the ionizing source, 
suggesting that quasars may have produced stronger helds 
than faint primordial galaxies. However, the latter conclu¬ 
sion is actually incorrect. Indeed, as we detail below, the de¬ 
pendence on the physical parameters is more involved than 
what equation (211 suggests. For example, increasing the 
luminosity of the source increases not only the number of 
photons efficient at yielding a larger B, but also the number 
of photons at the Hydrogen ionization threshold, leading to 
a bigger Stromgren sphere (though the dependency of Vs on 
Ltot is weak). Thus, in that case, the regions that get mag¬ 
netized are actually more distant from the source, and the 
geometrical dilution of photons must be properly taken into 
account in order to get the correct estimate of \B\. Hence 
a careful and detailed analysis of (121 is necessary, which is 
the purpose of the following section. 


4 APPLICATION AT COSMIC DAWN 

The analytical formula of the magnetic field ( |12[ ) we ob¬ 
tained is in principle applicable during Reionization in the 
vicinity of any ionizing source embedded in an arbitrar¬ 
ily inhomogeneous neutral medium. To gain further insight 
and obtain numerical estimates, we now apply this expres¬ 
sion to some simple models, considering specific sources and 
a mildly non-linear inhomogeneity outside their Stromgren 
sphere, as represented in figure 


4.1 Sources of ionizing photons 


We consider the three different types of sources that are 
believed to have driven Reionization: (clusters of) Popula¬ 
tion HI stars, primordial galaxies, and high redshift quasars 
(QSOs). For simplicity, each source is characterized by a 
power law luminosity, with normalisation Lq and spectral 
index a in a certain frequency range: 

L„ = Lo ^ ^ [^Oi ^i] (22) 


where we will call i/i the cut-off frequency. For quasars, 
we take a ~ —1.7, i^i = lOOt'o and a tota l luminosity of 
IO^^Lq which fixes Lq (e.g. Shang et al.|2011 1 . Values for the 
other sources were based on the curr ently available syntheti c 
spectra from the Yggdrasil modej^ (Zackrisson et al.|2011|. 


^ http://ttt.astro.su.se/~ez/ 



Figure 1. Representation of the modelled situation. A source 
forms a Stromgren sphere of radius r^. An inhomogeneity in the 
IGM is situated at a distance D from the source. All the graphs 
in this paper are plotted along the line of sight corresponding to 
S = ©max defined by equation l|34[|. 


which uses the Schaerer (20021 and Raiter et al. (20101 sin¬ 
gle stellar populations for population HI stars, differing by 
their initial mass functions (IMF). We used the PopHI.l 
model for Population HI clusters, a zero-metallicity pop¬ 
ulation with an extremely top-heavy IMF (50 — SOOM©, 
Salpeter slope), from which we obtained a ~ —0.3, — iuo 


sidered Population H stars with a Kroupa IMF and metal- 
licity Z G [0.0004,0.008], which yields a ~ — 2, = 4zzo 
and Lo = 3 X 10^® erg s“'^ Hz“^. In practice, we checked 
that our results do not depend on metallicity. From |Mar-| 
tini (20041, we assumed quasars to live about 100 Myrs. For 


Pop HI clusters and first galaxies, we considered constant 
star formation rates for 100 Myrs. 

For simplicity the Stromgren spheres of the sources are 
taken spherically symmetric. In reality, they are far from 
symmetric, exhibiting often a ‘butterfly’ shape in numerical 
simulations (e.g. Ciardi et al. 20011, so that = rs{9,ifi) 
in principle. In this case, the lower boundary of the inte¬ 
gral defining the column density in equation (181 varies spa¬ 
tially. In this first approach we neglect these angular varia¬ 
tions. The magnitude of the radius of Stromgren spheres 
depends on whether the recombination rate is sufficient 
to reach the steady state. The recombination time-scale is 
Lee = (asnC)”^, where n is the mean hydrogen density, 
as = 2.6 X 10“^^ cm® s“® is the case B recombination coeffi¬ 
cient at a gas temperature of lO’^ K and C{z) is the hydrogen 
clumping factor. This factor is still poorly constrained, and 
we adopt the fit C{z) = 27.466 exp(—0.1142 -|- 0.0013282®) 
obtained from simulations by Mellema et al. (20061. We note 
that the values of the clumping factor are somewhat sensitive 
to its definition (e.g. Finlator et al.|2012 l, but they remain 
of the same order of magnitude. In our redshift range, this 
yields free ~ 10 Myrs while our sources live 100 Myrs so, 
for simplicity, we considered the following expression of the 
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Stromgren radius 


f 3 iVion 


1/3 


(23) 


where the rate of ionizing photons emitted by the source 
is A^ion = Finally, outside Stromgren spheres, 

Thomson scattering is negligible since in neutral regions its 
cross section is by numerous orders of magnitude smaller 
than that of photoionization. 


z 

f I/O (kpc) 

liva (Fpc) 

hOvg (kpc) 

30 

0.0073 

0.47 

7.3 

15 

0.053 

3.4 

53 

10 

0.16 

11 

160 

6 

0.64 

41 

640 


Table 1. Orders of magnitude of mean free paths of photons of 
frequency i/ = 1/0,41/0 and IOpq at various redshifts during the 
EoR. 


4.2 Intergalactic medium 

For simplicity, we neglect the contribution of the first He¬ 
lium ionization and we consider a homogeneous Reioniza¬ 
tion scenario. We thus suppose that the ionization contrast 
Sx vanishes, meaning that the ionization fraction Xe is uni¬ 
form, so that Vxe = 0. In the rest of the paper, we will 
therefore focus only on the global term in equation ( | 12 [ ). 
The dynamical behaviour of baryons in the IGM is governed 
by 


dUe „ 2 dnm 

— =r,nH.-aan, = - — 


(24) 


where the photoionization rate at a distance r from the 
source is 

1 / O' j. —Tiy 

r;—Lue 




PCX 

J I'o 


hv 


'di/. 


(25) 


Outside the Stromgren sphere, since recombinations are neg¬ 
ligible, the variation time-scale of the densities tihi and rie 
is the ionization time-scale ti = T”'^, which depends on the 
distance from the source. These densities may be considered 
constant wherever ti/ts is smaller than one. The distance 
at which ti = ts is indicated by a vertical dashed line in 
hgures and To the right of this line, the assumption 
of constant nni and rie is perfectly safe. Between and this 
distance, their variation is not negligible in principle, but 
note that this line lies, in the vast majority of cases, very 
close to the Stromgren radius. Therefore the assumption of 
constant nni and Ue does not affect significantly the values 
of the magnetic field beyond this line. 

To model simply the inhomogeneous IGM in which the 
source is embedded, we consider an inhomogeneity next to 
the Stromgren sphere, centred at a position D from the 
source (cf. figure]^, with various prohles. The hrst simple 
profile consists in a Gaussian inhomogeneity 




riHi = n[ 1 + 5oe 2 <t 


(26) 


which is an overdensity for (5o > 0 and an underdensity 
for — 1 < 5o < 0. Such a Gaussian overdensity may result 
from gravitational instability, but inhomogeneities may also 
form through thermal instability, collecting the surrounding 
matter, in which case the overdensity is surrounded by an 
underdense region. Hence a second simple model of inhomo¬ 
geneity, which we will call hereafter a ‘Mexican hat’ profile 
(MH), composed of two imbricated Gaussians, a small width 
overdensity > 0 inside an extended width underdensity 


V < 0: 


whi = n l + SgC 


- 2 <t 2 




+ Sq e 


(27) 


Mass conservation f (whi — n) d^r = 0 then gives the con¬ 
straint (jo"(T+ -f (5(j"CT_ = 0. 


4.3 The global term 


In order to analyse in detail expression (16 1 and explicit its 


characteristics, we now consider an inhomogeneity with a 
Gaussian profile together with a source with a power law 
spectrum. Given the importance of the mean free path of 
photons in this model, we recall some typical orders of mag¬ 
nitude in table for frequencies typical of the ionizing pho¬ 
tons emitted by the sources participating to the EoR. 

4-3.1 The interaction term 


As we have seen in equation (16 1 , the strength of the mag¬ 


netic field depends on two criteria: the intensity of the in¬ 
teraction of photons at a given distance, quantified by Tint, 
and on the geometry of the situation, quantified by Egeom. 
Let us hrst focus on the former. For a power law spectrum 
( 221 , with the change of variable y = the interaction 


term reads 


= ^{r) 






7(5 


6-0 
T 3 

1 )- 7{5 


"o 


(28) 


where 7 is the lower incomplete gamma function. To lighten 
the expressions we have set ^(r) = | ('/o) ^ 0^0 

dropped the ‘global’ exponent from now on. 

Now dehne ri,g the distance from the source at which 
r,/o = 1, and r^^ such that — 1. Note that they roughly 
correspond to the mean free paths from the Stromgren 
sphere, r^. ~ +£ui, with equality for do ^ 1- Then, three 

regions may be dehned in the interaction term (cf left panel 
of figure namely: ‘very close’ to {rs < r r^g ^vi), 
‘close’ to {ri,g ^ r r^-h) and ‘far’ from (r„g <C r^^ ^ r) 
the Stromgren sphere. Given the smallness of (cf table 
[^, the ‘very close’ regime is not relevant for our purpose 
here as it lies within the width of the Stromgren radius. 


In the ‘close’ region, the interaction term (28 1 simplifies 


to 


F[ir = C(r) 


r(^) F(^) 


while ‘far’ from the Stromgren sphere 


« , « a —5 —T 

far \ / 1 ^0 \ { \ ^ 


F:nJ=C(r) 1 -^ 

Vl J \V0 J T, 


(29) 


( 30 ) 
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Figure 2. The properties of the magnetic field depend on whether the inhomogeneity is ‘close’ or ‘far’ from the Stromgren sphere 
(r <C or <S r, left panel), and on whether we are considering the regions ‘inside’ or ‘behind’ the inhomogeneity (r_ <C r <C r+ or 
r_|_ <g r, central panel). Thus, four different zones may be distinguished in which the magnetic field is the most significant (right panel). 
Hatching indicates regions where the strengths are weaker. 


In other words, the source has a certain ‘impact zone’ (the 
‘close’ region) inside which its ionizing photons are numer¬ 
ous enough to interact significantly. Outside this region, the 
strength of the field decreases exponentially. 


4-3.2 The geometric term 


Let us now focus on the second factor of expression (161. To 


see precisely what the requirement of a favourable geomet¬ 
rical condition consists of, let us consider the simple case of 


a Gaussian inhomogeneity (26l. In spherical coordinates as 


in figure the geometric term reads 


%/2(jn(5o 


Fi 


sin ft — 
a 


F2 


r — cos 6D 


-F2 


rs — cos 6D 


( 31 ) 


where 


and 


Fi{x) 



(32) 


^ 2 ( 3 ;) = e "2 - ^cose^ erf . (33) 

As far as the angular dependence is concerned, it is 
dominated by the Fi factor which is itself characterized by 
the two angles: 


Smax = arcsin 


Slim = arcsin 


:s 

3%/3 

2 



(34) 

(35) 


where Smax is the angle corresponding to the maximum 
of the function f’i(sinS^) and Sum is the angle for which 
the tangent at the inflexion point goes to zero. For angles 
larger than Sum, the strength of the interaction term is very 


strongly decreasing so that we will only consider lines-of- 
sight such that 9 £ [— Siim, Sum]- 

For the radial dependence, defining 

r±{6) = cos 6D F V2a, (36) 


three regions may be distinguished (cf. central panel of figure 
namely: ‘in front of’ (rs < r <C r-), ‘inside’ (r_ <C r ^ 
r+) and ‘behind’ (r+ <C r) the inhomogeneity. ‘In front of’ 
the inhomogeneity, the two F 2 terms essentially cancel each 
other out, so that B ~ 0 in this region. ‘Inside’ the geometric 
term simplifies to 


F 


inside 

geom 


nSoD^ 

- e 20-2 

2 a2 


sin 29 



while ‘behind’ the inhomogeneity 



FaJom'^ = ndoV^^sinS —e 
ra 


L eD)‘‘ 


(38) 


In other words, this geometrical factor decreases the 
strength of the field exponentially above a certain angle Sum- 
It hence constraints the magnetic field to be generated only 
in the vicinity and behind inhomogeneities, within a domain 
in the shape of a shadow. 


4-3.3 Magnetized regions 

The result is that the field has a rather simple spatial dis¬ 
tribution: It reaches its maximum strength in the vicinity of 
the inhomogeneity, roughly a.t 9 = Smax and r = D, and then 
decays with distance behind the inhomogeneity (cf figure]^. 
More precisely, to a given source corresponds a given de¬ 
limiting its ‘close’ and ‘far’ regions. Then the properties of 
the magnetic field generated around this source essentially 
depend on whether the surrounding inhomogeneities are in¬ 
side the ‘close’ or ‘far’ regions as depicted on the right panel 
of figure 
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Then we can see that if the inhomogeneity is close to 
the source (D < r^i), the strength reaches 




5-a\ 


-r 


c \ / T~\ \ s:f \ (o: —6)/3 

6 — a\ ( D — Ts + OQa \ 


, (39) 


while if it is far from it {D > r^j), it reaches 

_ / D — rs+5LcT\ 

□ inside ^ ni - A l'0'\ f I'o'] ) 

Bin. =tC(B)n5o(^l--j - 

(40) 

where Sq = Then after reaching these values, the 

field decays behind the inhomogeneity depending on the sit¬ 
uation according to 


□ behind -3 / ^ — Vs + 2Soa ^ ^ 

-^close ^ 


(41) 


S behind —3 

far OC r 


e^p(_r^n^) 

r-rs-l-25'<T 


(42) 


In these expressions we can see explicitly the role played by 
the amplitude of the inhomogeneity 5o- First of all, the sign 
of B is given by the sign of 5o. Also, we can see that the 
strengths are larger in underdensities than in overdensities. 
Indeed, for negative 5o, B -I- docr is smaller than for a 


positive do. Expressions (391 and (401 thus yield larger values 


for a negative do than a positive one. The same applies to 


expressions (411 and (421, showing that the fields will also 


be generated over larger distances in underdensities than in 
overdensities. 


4.4 Sensitivity to physical parameters 

Figure shows a typical result in the case of an underden¬ 
sity, for a primordial galaxy at 2 = 10. The superimposed 
lines are separating the different regions studied in the above 
section. Interesting field strengths are produced in regions 
that span transversally the domain between — Siim and 0iim, 
and extend radially over a few 14 ,^^ behind the inhomogene¬ 
ity. In addition, the dashed red arc represents the distance 
corresponding to ti = ts, below which tihi is overestimated 
as discussed in section [T^ The strength and scales reached 
by the generated field depend on many parameters: those 
characterizing the IGM [D, a or a+, So or and h), those 
characterizing the source (Lq, o and i/i), and the redshift 
2 . In light of the above analytical expressions, we now make 
those parameters vary one by one in order to evaluate their 
importance. 


4-4-^ Varying the properties of the inhomogeneity 

For simplicity all the configurations (source and inhomo¬ 
geneity) considered here possess an axial symmetry (cf fig¬ 
ure . Therefore there are no orthoradial gradients along 
the symmetry axis, and thus the magnetic field vanishes in 
the 0 = 0 direction. For this reason, all the graphs presented 
here are plotted along the line of sight 0 = 0ms.x- Note also 


- + 


- - 20.0 

- -20.5 

- - 21.0 

- -21.5 


- - 22.0 

-22.5 


Figure 3. Typical field generated by an underdensity at 2 = 10 
by a primordial galaxy. The thin black arc represents the edge of 
the Stromgren sphere. The thick black circle corresponds to the 
FWHM of the underdensity. The blue area corresponds to mag¬ 
netic fields pointing towards the reader, orange in the opposite 
direction. Strength is indicated in Gauss in logarithmic scale. Cor¬ 
responding to the analytical decomposition explicited in figure]^ 
the continuous red lines show the —6iini and ffuni directions, the 
dashed green lines correspond to the —6max and flniax directions, 
the blue arcs indicate r_ and and the continuous orange curve 
indicates = 4no here). The dashed red arc is the distance 

at which ti = ts (cf text). 




that all graphs start at r = r^. Obviously, the larger |(5o| 
the larger the gradients, and thus the higher the strengths 
generated in principle. On the other hand, large 5o imply 
strong absorption in the case of an overdensity. So for ex¬ 
ample behind a very dense region, the field vanishes but is 
important at its edges. In the following, we consider pock¬ 
ets of neutral gas, not virialized, in the vicinity of an early 
isolated ionizing source. Therefore we consider — 1 < 5o 5. 

In figure]^ we plot the field generated in an overdensity 
and an underdensity, differing only by the sign of 5o. We see 
that the strength is higher and generated on larger scales in 
the underdensity case, which is due to the fact that there 
are more photons at a given distance than in the overdensity 
case. This was already noticed previously from expressions 


(^ to 

In figure we study the influence of the position and 
profile of the inhomogeneity for a given source. For illustra¬ 
tive purposes we take a QSO at redshift 2 = 10. On the left 
of the figure, the profile of the inhomogeneity (a, do) is kept 
constant while its distance from the source D is increased. 
We can see that the closer the inhomogeneity, the stronger 
the strength, which is natural since photons are absorbed 
and diluted as they travel away from the source. Also, for 
QSOs, is huge (cf table so that we are always ‘close’ 
to the source and formulae and apply. We in¬ 

deed recover the power-law asymptotic trend with distance 
of equation ( |41[ ). In the central panel of figure D and 
So are kept constant as a varies. The narrower the density 
profile, the greater the strength, since gradients are more 
important for narrow profiles. Note that because the term 
we are considering in equation (121 is global, once a non 
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Figure 5. Magnetic field generated with various inhomogeneities, for a given source (here a quasar at z = 10). Left: Gaussian overdensity 
of fixed profile, varying distance. Middle: Gaussian overdensity at fixed position, varying its width cr. Right: Mexican hat profile at fixed 
position, varying its width (T+. The lines corresponding to distances at which ti = tg (cf. section [4.2[ | are not visible here because for 
quasars the Stromgren sphere is so large that photons are very diluted as they reach rs, and ti becomes much larger than tg already at 
the edge of the Stromgren sphere. 



r(kpc) 


Figure 4. Magnetic field generated by a primordial galaxy, with 
two Gaussian inhomogeneities of same |<5o|, but with opposite 
signs (blue: positive; black: negative). We see that the mecha¬ 
nism is more efficient in underdense regions. Vertical dashed line 
indicates the distance at which ti = ts- 


radial gradient is formed, it generates magnetic fields along 
the entire line of sight. This is why the narrow profiles gen¬ 
erate fields of stronger strengths on larger distances. Their 
angular extent about the symmetry axis is smaller though. 
The right panel of has to be compared with its central 
panel: They correspond to the same study but with a MH 
profile instead of a Gaussian overdensity profile. We can see 
that the strengths are larger. This comes from the fact that 
photons in underdensities are less absorbed so that the flux 
of ionizing photons is more important within the overden¬ 
sity than in the simply Gaussian overdensity case. This is all 
the more true than the width of the underdensity is impor¬ 
tant, which is why contrary to the gaussian case, the wider 
the profile, the larger the strength for an MH profile. Also, 
since the field has opposite directions in underdense and 
overdense regions, the contributions from the underdensity 
and the overdensity composing the MH profile cancel out. 
This is why for such a profile, at some distance, the field 
changes sign. In the right panel of figure we can see that 
this distance depends on the width of the MH. 

In both figures and 1^ when the inhomogeneity is very 
close to the Stromgren sphere, the strength reaches high val¬ 
ues at the outer edge of the ionized region. This is because in 
our model we impose spherically symmetric ionized regions 
which intercepts the Gaussian profile of the inhomogeneity 
thus inducing non radial gradients and a magnetic field. In 
more realistic situations, such border effects will indeed take 


place as the ionization front will be distorted at the contact 
of an inhomogeneity. A more precise study of such effects is 
left for future work. 

Finally, expression ( |18[ l dictates the topology of the 
magnetic field. Field lines given by B x dZ = 0 in spheri¬ 
cal coordinates with dl = dr r + rdO 9 -\- r sin Odcf) (j> satisfy 


dr = Q and 


de 

dtp 


dip J riHidr 

de f riHidr ' 


(43) 


Since dr = 0, field lines remain on spheres centred on the 
source. In the case of a (^-independent configuration, as in 
all the examples considered here, dd = 0 so the lines are 
circles around the axis of symmetry of the system. This is 
why in all the density plots of this paper the fields point 
perpendicularly to the plane of the plot. In a more general 
case the lines remain loops at given radii since dr = Q but 


with a more complicated shape given by (431. Though, once 


formed, these fields will be processed by the velocity field of 
the IGM, their peculiar initial spatial configuration is inter¬ 
esting as it may help discriminate magnetic fields generated 
through this process from those generated by different mech- 


4-4.2 Varying the properties of the source 

In figurej^we choose a certain MH inhomogeneity, at 2 = 10, 
and vary the spectral index a and cut-off frequency Vi of the 
source, here a primordial galaxy for illustration. 

On the left of the figure, we make the cut-off frequency 
vary. We observe that above a certain value, the strength 
of the generated magnetic field has reached a plateau. For 
instance in the case shown, the field strength generated 
with = lOr'o is almost identical to that generated with 
zzi = 3000r'o- This is because iiovg larger than 50 kpc at 
2 = 10, so photons of frequency above IOzzq do not inter¬ 
act significantly in the neighbourhood of the inhomogene¬ 
ity. Also, far behind it they are too diluted to generate a 
significant magnetic field. This result underlines the follow¬ 
ing compromise: This mechanism requires photons of high 
enough energy to photoionize deep in the IGM to generate 
magnetic fields on large distances, but not too high either 
because of dilution (cf. equations ( |41[ ) and ( |42[ |). 

On the right of figure we make the spectral index a 
vary. Since the total luminosity is fixed, harder spectra cor¬ 
respond to sources with fewer low energy photons, which is 
why harder spectra generate weaker fields near the source. 
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Figure 6. Magnetic field generated with various source properties, with a given Mexican hat inhomogeneity at a given position, at 
z = 10. Left: For a fixed power-law spectrum, we vary the cut-off frequency ui. Right: For a fixed total luminosity, we vary the spectral 
index a. Vertical dashed lines show the distances at which ti = tg. 


but stronger fields far from the source. Also, naturally, the 
harder the spectrum the further in the IGM the field is 
generated. However, we can see in this example that the 
strengths do not depend too strongly on the spectral index, 
only up to factors of two. This can be seen in the analytical 
expressions of the previous section like equation (281, where 
(a — 5)/3 and {a — 6)/3, for the relevant values of a, do not 
vary much with a and so B does not either. 

Finally, we consider various redshifts. Depending on the 
epoch, the type of ionizing sources and the properties of the 
IGM are different. For illustration, we decompose the EoR 
in three stages: PopIII star clusters for z £ [30,20], primor¬ 
dial galaxies for 2 : £ [20,10] and quasars for z G [10,6]. 
Figure]^ shows examples of field configurations obtained at 
z = 20 and a = 10, and table is a summary of the typical 
strengths obtained after varying the different parameters of 
the model. They allow us to identify the following trends. In 
the beginning of EoR, Population III clusters generate rel¬ 
atively high strengths (~ 10~^® G) but on relatively small 
scales (hundreds of parsecs), while on the contrary, at the 
end of EoR, quasars generate low strengths (~ 10“^^ G) 
but on important scales (hundreds of kpc). First galaxies 
are somewhat in-between and generate ~ 10“^^ G on tens 
of kpc in the middle of the EoR. This may be interpreted 
as follows. Population III clusters are not very energetic so 
they form small Stromgren spheres and thus photons from 
the source are not too diluted as they reach the IGM. How¬ 
ever, since they appear early in the EoR, n is large and thus 
the mean free path of photons is small. Therefore many pho¬ 
tons photoionize and generate strong magnetic fields near 
the Stromgren sphere. On the contrary, quasars are very en¬ 
ergetic and generate huge Stromgren spheres, so that pho¬ 
tons from the source are very diluted as they reach the IGM. 
In addition since quasars appear late in EoR, the IGM is not 
very dense and the mean free path of energetic photons is 
huge. Consequently the ionizing flux is relatively low and the 
generated magnetic field is weak, but extends on important 
scales. 


5 DISCUSSION 

The cases detailed above are idealized. In reality, as we have 
mentioned, Stromgren spheres are deformed, there is a full 
distribution of asymmetric inhomogeneities outside the Hll 
regions, and ionizing sources may not be isolated. Taking 
all this into account would allow us in principle to compute 


Source 

Redshift 

Log |B| 
(Gauss) 

Distance from the 
ionization front (kpc) 

Pop III 

30 

-19 

0.3 



-21 

1 


20 

-19 

0.5 



-21 

1 

Primordial 

20 

-20 

10 

galaxy 


-22 

15 


10 

-21 

30 



-22 

100 

Quasar 

10 

-21 

300 



-22 

1000 


6 

-22 

500 



-23 

1500 


Table 2. Third column shows the typical strengths generated at 
the distances away from the Stromgren sphere shown in the fourth 
column. These strengths are typical values obtained by varying the 
different properties of the inhomogeneities. 


the full magnetic power spectrum generated by photoion¬ 
ization. Furthermore, the IGM is dynamical and the weak 
turbulence associated with the mildly non linear regime of 
structure formation on large scales will process the initial 
power spectrum of the magnetic field. A study of all these 
effects is beyond the scope of this paper and should be best 
addressed with dedicated numerical simulations. 

As we have seen, an interesting feature of this mecha¬ 
nism is that it generates magnetic fields outside Stromgren 
spheres on distances of the order of a few , or even larger 
in underdense regions. By comparing these distances to the 
half typical separation between Stromgren spheres, we may 
get an idea of the fraction of the Universe that may be mag¬ 
netized by this mechanism. Let us first do so for Popula¬ 
tion HI star clusters at 2 = 20. Considering they formed by 
molecular cooling in gas overdensities of mass 10®M© (c.f. 
Barkana fc Loeb(2000 1, their half physical mean separation 


is roughly 10 kpc, estimated from the abundance of their 
~ 3a parent haloes (Mo & White 20021. In addition, the 
radius of their Stromgren spheres is of order a fraction of 
kpc, and £ 4^0 is about one kpc large. Thus, it seems that 
these sources magnetize essentially on the outskirts of their 
Stromgren spheres, leaving an important fraction of the IGM 
unmagnetized. Consider now faint dwarf galaxies at 2 = 15, 
and quasars of IO^^Lq at 2 = 10. Faint dwarf galaxies are 
the major candidate sources responsible for reionization (e.g. 
[Bouwens et al.|[2012[ [Wise et al.||2014| . Assuming they are 
hosted in 10“ Mq haloes, their half physical mean separation 
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PopIII at z = 20 



|B(G)| 


|B(G)| 


First galaxy at z = 10 

hhi (cm"-'*) 




Quasar at z = 10 


nHi(cm ■'') 




kpc 


Figure 7. Examples of the field generated with Mexican hat inhomogeneities for each type of sources considered here. Left column: 
Left y-axis and blue curves show the absolute value of the magnetic field strength. Vertical dashed lines indicate the distances at which 
ti = ts. Right y-axis and dotted lines show the inhomogeneity profiles. Right column: Black line is the Stromgren sphere. Inner and 
outer dashed black circles correspond to (j+ and cr_ of the inhomogeneity, respectively. Blue areas correspond to magnetic field pointing 
towards the reader, orange in the opposite direction. Strengths are in Gauss in logarithmic scale. 
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is a couple of tens of kpc. Their Stromgren radii are of the 
order of a few tens of kpc, and £ 4^0 is a couple of kpc large. 
Luminous quasars, on the other hand, are extremely rare at 
2 ; = 10, which seems confirmed by the recent Planck results 
implying that the redshift of reionization, Zre = 8.8, is lower 


than previously thought (Planck Collaboration 20151. For 


illustration, considering they are hosted by 5a haloes, their 
half physical mean separation is of a couple of Mpc. Their 
Stromgren spheres have radii of the order of the Mpc, and 
these sources have photons of large enough mean free paths 
to magnetize in between these spheres. Therefore, these or¬ 
ders of magnitude are more favourable in the case of pri¬ 
mordial galaxies and quasars than in the case of Population 
III star clusters, and suggest that these sources may have 
participated through this mechanism to the weak magne¬ 
tization of an important fraction, if not the whole, of the 
IGM. Note that even the case of Population III star clusters 
is interesting, as these sources thus premagnetize the envi¬ 
ronment in which the next generation of stars and galaxies 
forms. 

In this study we have focused on the photoionization 
term in the induction equation. An exact comparison of this 
term with the Biermann term is beyond the scope of this pa¬ 


per. We note that Doi & Susa (20111 examined their relative 


importance in numerical simulations of the neighbourhood 
of an ionizing super-massive star. In their study, they focus 
on the situation across an ionization front, in which a self- 
shielded, neutral, 5 ~ 10^ — 10® over-density defines very 
sharp and strong gradients in the temperature and elec¬ 
tronic density fields. Such a situation could indeed occur 
within Stromgren spheres of the very first luminous sources. 
Under those conditions, they concluded that the Biermann 
battery dominates by one order of magnitude. In our case, 
we considered mild, neutral over-densities way outside the 
Stromgren regions of stronger, long-lived ionizing sources. In 
such contexts, the Biermann battery may not be effective, be 
it for purely geometrical reasons. To see that, consider again 


equation (111, in which we have separated the global pho- 


toionizing contribution from the local contribution. It is the 
global contribution that creates magnetic fields in regions 
where the Biermann battery does not operate. Comparison 
of the two mechanisms is relevant only in regions where they 
coexist. For that purpose, it is convenient to note that equa¬ 


tion (111 may be rewritten as 


as = 5^ 

e ni 


X [-Vpe +Ppe\ - X Ppe. 

L J eue 


(44) 


On the right hand side, the last part contains the full global 
term examined in this paper (cf. equation |16[ ), and a con¬ 
tribution to the local term. The remaining contribution to 
the local term is contained in the Ppe present in the square 
brackets. In regions where the Biermann battery is efficient, 
i.e. where pressure and density gradients are not aligned, 
it is enough to compare Ppe to the pressure gradient. More 
precisely, since Ppe is always radial, their comparison makes 
sense only in cases where —Vpe is radial too. Then, consid¬ 
ering a perfect gas equation of state Pe = n^khTe, their ratio 
yields typically 


Ppe 




13 


Te 

10 K 


10 kpc 


(45) 


where Lg is the typical scale of temperature gradients. For 


illustration, we have taken ppe = 3 x 10“"^^ erg cm“"^, which 
is obtained at a rather remote distance of 2rs ~ 100 kpc 
from a primordial galaxy at 2 = 10. Hence, in this case, the 
photoionization term dominates, even at a distance of 2rs at 
which photons are very diluted. Also note that, because the 
Biermann battery contribution is independent of the prop¬ 
erties of the source, the ratio above scales linearly with the 
luminosity of the ionizing source. 

The magnetogenesis mechanism examined here, based 
on the simple physics of momentum transfer from ionizing 
photons to photoelectrons, is therefore likely to have gen¬ 
erated seed magnetic fields on cosmological scales during 
Reionization. Although weak and very remote, the strengths 
of the seeds produced, together with their specific spatial 
configuration, could actually be revealed directly through 
the recently proposed probe of magnetic fields in the EoR 
detailed by Venumadhav et al. (20141, although large coher¬ 
ence lengths of the magnetic fields might be mandatory. 
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APPENDIX A: MOMENTUM TRANSFER 
TERM 

A fraction fmt of the momentum of the incident photon is 
transferred to the freed electron during photoionization: 

meF=r, (Al) 

c 


where /mt 

TgM: 


is frequency dependent (Sommerfeld & Schur 


fmtO = 


8 u — Vo 

5 V 


(A2) 


Note that this fraction may be larger than one, in which 
case the ions recoil. 


Now, by definition dtfe\^d^vd^rdt is equal to the number 
of photoelectrons of speed v in direction v generated at a 
position r at a time t. Since we consider Hydrogen only, each 
photoionization produces only one electron. This number is 
thus equal to the number of photoionizations due to photons 
of frequency v in direction k — v where v satisfies (All. 
Finally, since the photoionization rate density is the product 
of the number density of incident photons, the velocity of 
incident photons, the number density of target Hydrogen 
atoms and the cross section, we get 


dtfe\^d^vd^fdt= ccr^nHid^rdt (A3) 

where the number density of incident photons of frequency 
V with direction A: at r at time t is 

= (A4) 

by definition of the monochromatic specific intensity. There¬ 

fore we model the source term microscopically in equation 
0by 

dtfel^d^v = ^^j^nmdvdid, (A5) 

so that macroscopically we get from expression from 

equation (0. 


APPENDIX B: NOTATIONS 


In the multifluid description, macroscopic physical quanti¬ 
ties are dehned as 


nc, = f fad^v 

Va = J V fad?V 

Pa = TTlaf {Va “ (^14 “ fad^V 

J OL — 


(Bl) 


respectively the number density, the velocity, the pressure 
tensor and the current density of species a. This description 
is then reduced to a single fluid with 


Pm — 


Pq — 


Sq. T 1 

J - j cx . 


r> _ v jdCM 

Jr — Zjq.J Q, 

Pf" = / (U - f) (U - f) fad^v 


(B2) 


respectively the mass density, the charge density, the centre- 
of-mass velocity, the current density and the total centre-of- 
mass pressure tensor in the one-fluid. Then, taking the first 
moment of 0 weighted by qa , and summing over all species, 
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yields the generalized Ohm’s law Q, with the additional 
notations: 


P = T.o 


-V • p: 


.CM 


C = T^cQc J V dtfd^d^v. 


(B3) 
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